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We analyze the relationship of longitudinal and transverse increment statistics measured in 
isotropic small-scale turbulence. This is done by means of the theory of Markov processes lead- 
ing to a phenomenological Fokker-Planck equation for the two increments from which a generalized 
Karman equation is derived. We discuss in detail the analysis and show that the estimated equation 
can describe the statistics of the turbulent cascade. A remarkably result is that the main differences 
between longitudinal and transverse increments can be explained by a simple rescaling-symmetry, 
O |' namely the cascade speed of the transverse increments is 1.5 times faster than that of the longitudi- 

D . nal increments. Small differences can be found in the skewness and in a higher order intermittency 

. term. The rescaling symmetry is compatible with the Kolmogorov constants and the Karman equa- 

tion and give new insight into the use of extended self similarity (ESS) for transverse increments. 
Based on the results we propose an extended self similarity for the transverse increments (ESST). 
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I. INTRODUCTION 



PACS numbers: 05.10.Gg, 47.27.-i, 47.27.Eq, 47.27.Gs, 47.27.Vf 

i 

CO 

. $h '. Small scale turbulence is not yet full understood. A complete theory based on the Navier-Stokes equation has not 
been achieved yet, thus our present understanding relies for the most part on phenomenological and experimental 
approaches. It is assumed that turbulence forms an universal state which exhibits stationarity, isotropic and homo- 
geneity in a statistical sense j2^,|44|. In general, turbulence are driven on large scale, i.e. energy is injected into large 
scale motions and is dissipated on small scales resulting in a net flux of energy from large to small scales j5||. The 
energy flux results from the inherent instability and the subsequent breakup of vortices into smaller ones. 

For local isotropic turbulence the main challenge is to understand spatial correlation. Usually the turbulent field 
U(x, t) is characterized by increments u r = e • [U(x + r, t) — U(x, t)], where e denotes a unit vector in a certain 
direction, x denotes a reference point and r a displacement vector. The increments are taken as stochastic variables 
in dependence of the scale variable r = |r| and by varying r correlation on different scales can be studied with 
0^ ■ the increments. In the following, we assign u r to longitudinal increments, that means e is parallel to r and v r to 
transversal increments, for which e is orthogonal to r. For specific length we write Ui and Vi in a short way. 
The central challenge in turbulence is to explain the statistics of exceptional frequent occurrence of large velocity 
increments on small scales r, which can not be understood with normal statistics. This is the so called intermittency 
CO ( problem. 

The work of Kolmogorov [33, is still the foundation of the small-scale turbulence theory. The starting hypothesis 
is that the possible symmetries of the Navier-Stokes equation are recovered in a statistical sense for high Reynolds 
numbers. These symmetries are homogeneity, i.e. the statistics of the increments is independent of the reference 
point x, isotropy, i.e. the statistics does not change under rotation of the frame of reference, and stationarity, i.e. the 
statistics is not time dependent. A further hypothesis is scale invariance, i.e. loosely spoken the structures of different 
sizes looks similar. Kolmogorov has considered the statistics in terms of structure functions (moments of the velocity 
increments), for which scale-invariance reads as (it") = (r/r{) (it™) implying (it™) cx r^™ with two different distances 
r and n. 

Using this hypothesis, Kolmogorov has furthermore assumed that for high Reynolds number the statistics of velocity 
increments depends only on the energy dissipation e and the scale r and ended up by using dimension arguments 
with (u™) = C„e™/ 3 r"/ 3 for r\ -C r -C L. The constants C n denote the Kolmogorov constants, r] is the scale where the 
dissipation takes place and L, the integral length scale, is the largest scale of the flow. However there are correction 
due to the fluctuating dissipation energy leading to deviations from the exponent £ n = n/3. This typical property of 
turbulence is an other aspect of the above mentioned intermittency. According to the refined self similarity hypothesis 
(RSH) of Kolmogorov, which takes into account the fluctuating energy dissipation (e r ) averaged over a volume V with 
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the extension r (e r = J v edx), the scale-invariance is given by 

«) =c„(e; i/3 )r"/ 3 cxr«". (1) 

The first model for the exponents is Kolmogorov's log-normal model, which results in = rt/3 — /ttn(n — 3)/18 
|39j |. If instead of the structure functions the probability density functions (pdf) p(u r ) are considert, this nontrivial 
scaling behavior correspond to a change of the form of the pdfs with r. To explain intermittency, i.e. the non-trivial 
scaling behaviour, is still one of the prominent challenges in turbulence. 

At first, the examinations were only concentrated on the longitudinal increments. There was hope that the exponent 
£ n can describe the self-similarity of every quantity of the velocity field. But recently, a lot of afford was undertaken 
to consider also the transverse increments and they seem to have essential different properties Q, |E El IE El El El 

E3, EE EE EE EE 123, m, H3, HH H3, [S3, m, 113, E3, EE E3, EE H3, HE E3, IH, E3l Thi s implies that the models as 

well as the considerations concerning only the longitudinal increments are too specific and a better understanding of 
the turbulence has to include the transverse increments. 

Longitudinal and transverse increments belong to different geometric / kinematic structures in the flow. The 
longitudinal increments can be associated with strain-like structures, the transverse increments with vorticity-like 
structures. Thus it is natural to modify RSH of equation and to use local averaged enstrophy (squared vorticity) 
rather than local averaged dissipation for the transverse increments (refined similarity hypothesis for transverse 
increments, RSHT, see [lUll5l p. 

«) = C t , n (fi^ 3 )r"/ 3 ocr^". (2) 

If intermittency results from the fluctuating dissipation and enstrophy, then the deviations from the Kolmogorov 41 
theory [37j can be investigated by the scaling of (e"^ 3 ) and (J7^ 3 ), respectively. For infinite high Reynolds number 
the scaling of averaged dissipation and enstrophy should be the same Eg, but in many experiments and simulation 
one observe differences in both exponents for finite Reynolds number fill IbH ITsl El EJ EE El EE El El. At 
least four arguments are used to explain this observations: 1) anisotropy, 2) effect of Reynolds number, 3) influence 
of boundary conditions, 4) intermittency. Anisotropy typically exists in every flow on large scales and it can influence 
the exponent stronger than intermittency itself 58, 69]. Decreasing the Reynolds number effects the scaling exponents 
in such a way that the differences between them increase EJ, l58t [69 j . 

The structure functions (u™) and (u") as well as the scaling properties involved can not define unambiguously the 
turbulent field; for instance flows with different structures may have the same scaling properties [64^. The same is 
true for the probability distributions p(u r ) and p(v r ), which are in essential equivalent to the structure functions. 
The reason is that these quantities are just two-point statistics. Definitive more general and detailed are the multi- 
point (or multi-scale) distributions p(ui, v\, 7*1; u n ,v n , r n ) (or multi-scale structure functions (J\ i u™ 1 Y[j v J l3 )) f° r 
different scales r,. These probabilities also consider the joint statistics of longitudinal and transverse increments and 
thus enables also to describe the interaction between them. Furthermore they describe the simultaneous occurrence 
of increments Ui, Vi on several length scales 

In this paper, we focus on an approach to characterize multi-scale statistics. It has been shown that it is possible 
to get access to the joint probability distribution p(u±, n; W2, Ti\ . . . ; u n , r„) via a Fokker-Planck equation estimated 
directly from measured data [20I I2H E3 • This has attracted interest and was applied to different problems of the 
complexity of turbulence like energy dissipation EE EE EE > universality of turbulence E^ , the theoretical derivation 
of the Fokker-Planck-equation from the Navier-Stokes equation in the limit of high Reynolds number [17j and the 
analysis of stochastic time series EE EE E2- The characteristic of this method is that it is based on pure data 
analysis, i.e. it is a parameter-free method, and that the few underlying assumptions are verifiable. Thus no special 
model-ideas are interwoven with the procedure. 

In this article we extend this approach to analyze the joint statistics of longitudinal and transverse increments. 
A first result concerning the different cascade speed of longitudinal and transverse cascade has been published in 
a preceding letter 60]. The aim of this article is to present the used method and its extension to longitudinal 
and transverse increments in detail and to discuss the similarities and differences of longitudinal and transverse 
increments. We start with a short description of the concepts of Markov processes and define the notation. In section 
IIIII wc describe the experiment and characterize the data with two-point statistics (i.e. with one scale) in section llVl 
Next, we analyze longitudinal and transverse increments separately by means of Markov-processes with respect to 
multi-points or multi-scale statistics. Then we go over to a combined analysis of both increments and discuss a new 
symmetry which connect both increments. An interpretation in section Ivl and a conclusion will finish the paper. 
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II. THE MATHEMATICS OF MARKOV PROCESSES 



The basis of our argumentation and of the used analysis is the theory of Markov and diffusion processes. Therefore 
we will give a concise description of them, which also serve as a guideline for the analysis of the turbulent signal. The 
foundations are known since Kolmogorov |36j |. but a detailed presentation can also be found in standard textbooks 
such as [HIHEtI 

We restrict ourself to the case of a two dimensional stochastic variable, denoted by the stochastic vector 

^ = ( Z ) • (3) 

Usually, a stochastic process is formulated in the time t but one can also imagine different independent scalars; for 
our purpose to characterize the turbulent cascade the scale r is the independent variable. 

The stochastic process underlying the evolution of u r in the scale r is Markovian, if the conditional pdf 
p (ui, ri | U2, Ti \ ■ ■ • ; u n , r n , ) with r\ < T2 < ■ • • < r n fulfills the relation 

p(ui,n | U2,r 2 ;u3,r 3 ; . . . ;u„,r„ ) = p(m,n | u 2 ,r 2 ) , (4) 

where p(ui,ri | u 2 , r 2 , U3, r 3 ; . . . ; u„, r n ) denotes the probability for finding certain values for Ui at some scale n, 
provided that the values of at the larger scales r^ > ri, i > 1, are known. The condition Q states that the 
increment distribution of U\ on r-y only depends on the increment value on one larger scale r 2 and that further scales 
do not give more information. Markov properties imply a remarkable property. The knowledge of the conditional pdf 
p ( u r , r I uo, rg ) (with r < ro) is sufficient to determine any n-point pdf: 

p(ui,ri, u 2 , r 2 , . . . , u„, r n ) = p (ui, n | U2, ri ) X p (112, vi \ 113, r% ) X ... 

Xf>(u„_i,r„_i |u„,r„) xp(u„,r„), (5) 

i.e. the single conditioned pdfs contain the entire information about the stochastic process 

For Markov processes the evolution of the conditional pdf in the scale r can be described by the Kramers-Moyal 
expansion, a partial differential equation for p ( u r , r | Uo, ro ) in the variables u r and r. According to Pawula's theorem, 
this expansion truncates after the second term if the fourth order expansion coefficient vanishes [ST!- In this case, the 
Kramers-Moyal expansion reduces to the Fokker-Planck equation (or Kolmogorov equation) 

-r-^p{u,r\u 0l r a ) = ( D^^u.r) p(u,r | u ,r ) ) (6) 



=1 
2 



d 2 



( D\f(u, r) p( u,r\u a ,r ) 



Note that we have multiplied both sides of the usually used Fokker-Planck equation with r; on the right hand side 
we have incorporated this factor in the definition of the Kramers-Moyal coefficients (JHJ . The Fokker-Planck equation 
describes not only the evolution of the conditional probability p( u r , r \ u ro , ro ) but also the evolution of the pdf 
p(u r , r), as can bee seen by integrating JHJ over uq. 

The drift vector and the diffusion matrix D^ 2 ) of ®) are defined via the limit 

£> 4 (1) (u,r) = lim M^\u,r,Ar), 

D^\n,r) = lim Mg\u,r,Ar), (7) 

where the coefficients M( fc ) are given by 

M^Vu, r, Ar) = — ( (u'Ar - Ar) - uAr) ) I u,r ) , 
Ar 

Mg\u,r,Ar) = -1- ( ( u J( r - Ar) - m(r)) x 

x (uj(r- Ar) -Uj(r))\u,r ). (8) 

The coefficients M^^ are conditional expectation values (-|u) on the veloctiy increment u and can easily be determined 
from experimental data. One can find estimates for the by extrapolating the measured conditional moments 
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]V[( fe ) in dependence of Ar towards Ar = [54[ . An other possibility would be to approximate the drift and diffusion 
coefficients by the coefficients M( fc ) for one finite Ar. The corrections of higher order in Ar can then be taken 
into consideration by correction terms. The approved coefficients can be used to recalculate the corrections, and 
if one perform this procedure recursively, the result converge towards the drift and diffusion coefficients |53j| . Both 
approaches give similar results, but the approximation of [53| has some pitfalls, if it does not examined the limes, 
which is an essential point of this analysis, see also (2^ • 

A solution p(u, I) of the Fokker-Planck equation can be derived from the Chapman-Kolmogorov equation p(u, I) = 
J p(u, 1\uq, Io)p(uq,Iq)(1uq. Here we use the logarithmic scale / := ln(i/r) to transform the Fokker-Planck equation 
into the usual form dp/dl = . . . instead of —rdp/dr—. . . . For sufficient small Al = I — Iq, such that are constant 
in Al, the conditioned probability can be approximated by 

p(u,l + Al\u',l) « 1 =x (9) 



47iVdetD( 2 )(u',l)AZ 



x exp 




(d^V))" 1 ] (u 3 ^ Df Al) (u k u' k D^Al) 

ijk 



For larger distances I — Iq, a solution can be obtained by reiterating the Chapman-Kolmogrov equation with this 
approximation: 

p{u,l) = J dujv-i ... J du p(u,r|u W -i,ijv-i) • ■ -Jj(ui, h\u , lo)p(uo, l ). (10) 

Below we use this path-integral approximation to solve the Fokker-Planck equation numerically. 

The error analysis for a diffusion process in two variables is more complicated than for one dimension. We perform 
the error analysis according to |6lj; here we just sketch the procedure and give the resulting equations for the error 
estimations. An error analysis comprises a quantity to estimate and a stochastic quantity which leads to an uncertainty. 
The two quantitates to estimate are the drift and diffusion coefficients and it is assumed that the stochastic quantity 
is given solely due to the stochastic nature of the process. Because the randomness is determined by D^ 2 \ the error 
of L>W an d £>( 2 ) can be expressed by itself. The following results are valid for small Ar. The error of the drift 
coefficient is given by 



A^ (1) = ^=V^\ (11) 

where N is the number of samples used to estimate D^ 1 '. To calculate the error for the diffusion matrix, one has 
to transform the diffusion coefficient by a suitable orthogonal transformation matrix B in the diagonal form. In the 
diagonal system one has 



D« = BD«, D< 2 > = BD^B T (12) 



— (2) 

The error of is then given by 



AD^ = D$a(m = D^/Jd^,N), (13) 



which leads to the error of the initial diffusion matrix by the inverse transformation AD< 2 ) = B T AD( 2 'B. The 
function a(m,N) is defined by 

<0+<T 

J P (CR = 0.68, (14) 

where 

N 

N 



~ 1 N 



»=i 



is an estimator for the average value of a x 2 -distributed stochastic variable £ = (to + T) 2 with a normal distributed 
stochastic variable T obeying (r) = and (r 2 ) = 1. Thus a(ni, N) defines the 32%-confidence interval of the estimator 
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In the above presentation the central quantity was the pdf but very often one is interested in the structure functions 
(u m v n ) — J u m v n p(u,v,r)dudv. From the Fokker-Planck equation a hierarchical equation for the structure functions 
can be derived by using (u m v n f(u,v)) = J u m v n f(u,v)p(u,v 1 r)dudv: 

-JL/ U m v n ) = (16) 
or 

+ m(u m - 1 v n D^ ) (u,v,r)) + n(u m v n - 1 D { 2 ) (u,v,r)) 

+ mn(tt m -V- 1 D$(u,t; > r)). 
This equation allows a direct comparison of the Fokker-Planck equation with structure functions. 



III. THE EXPERIMENTAL SETUP 



For the subsequent analysis we use two data sets measured in the central region of a wake behind a cylinder. The 
Reynolds numbers are R\ = 180 and R\ = 550, respectively. If nothing else is said, the data set with R\ = 180 is 
used. The high Reynolds number data set is used for comparison. 

For the first data set the local velocity is measured in a wake 60 diameters behind a cylinder with cross section 
D = 20 mm. The Reynolds number is 13200 with a Taylor based Reynolds number of R\ =180. We have measured 
the velocity component U in the mean flow direction, the V component orthogonal to the cylinder axis with an X-wire 
(Dantec's frame 90N10 with Dantec's 55P51 X-wire). We collect 1.25T0 8 samples with a sampling frequency of 25 kHz 
using a 16bit A/D converter; high frequency electronic noise was suppressed with a low-pass filter to prevent aliasing. 
We use Taylor's hypothesis of frozen turbulence to convert time lags into spatial displacements. With the sampling 
frequency of 25 kHz and a mean velocity of 9.1 m/s, the spatial resolution of the measurement is not better than 0.36 
mm. For the integral length scale L — R(r)dr, where tq represent the first zero-crossing of the autocorrelation 
R(r) function, we obtain a value of L =137 mm for the stream- wise component U and L t =125 mm for the component 
V. Taylor's microscale A is calculated using the method proposed by J5J- For isotropic turbulence, A can be written 
as 

where a is the standard deviation of the turbulent fluctuations. The limit has been calculated by fitting a second-order 
polynomial to the data, resulting in A =4.8 mm for the stream- wise component and At =3.0 mm for the V^-component. 
The dissipation scale ij can not be resolved, because it is smaller than the length of the hot-wire and smaller than the 
spatial resolution calculated with the finite sampling frequency and the Taylor- hypothesis. 

The high Reynolds number flow is measured 40 diameter behind a cylinder with a cross section of 5 cm. The 
Reynolds number is Re= 84600, the Taylor based Reynolds number is R\ — 550. The [/-component is measured in 
direction of the mean flow, the V^-component in direction of the cylinder axis. 1.25 • 10 7 data points are measured 
with a frequency of 200 kHz. The facilities are the same as in the first measurement. Typical length are L — 15.5 cm 
and A = 3.87 mm 



IV. EXPERIMENTAL RESULTS 



A. Two point statistics 



First of all, we characterize our data by two-point statistics (structure functions, pdf of increments etc.). In this 
way, we have a reference to results presented frequently in the literature and we get a first impression of the differences 
between longitudinal and transverse increments. 

A central pre-assumption for many considerations of small scale turbulence is isotropy. To study isotropy we use 
the Karman equation (isotropy relation) because it is the simplest exact relation for structure functions which rely 
on isotropy. The Karman equation connects the second order longitudinal and transversal structure functions: 

v d 

( v r) — ( u r) + oo - ( u r) nrst Karman equation. (18) 
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FIG. 1: Second order transverse structure function estimated on the one hand directly from data (squares) and on the other 
hand from the longitudinal structure function via the Karman equation (straight line); i.e. the left and right hand side of 
equation l|18|l are plotted. The relative deviations are plotted by a dashed line. For smaller scales r < L the Karman equation 
is well fulfilled within 5% and is isotropic in this sense. On larger scales the data are clear anisotropic. 

Fig. n]shows the left and right hand side of the Karman equation, i.e. the second transverse structure function and the 
hypothetic transverse structure function for isotropic flows calculated from the longitudinal one. Taking the validity 
of the Karman equation as indication of isotropy, we see that for small distances r < L the isotropic relation is well 
fulfilled, the relative error is within 5%. For large distances r > L isotropy is violated (large scale anisotropy), the 
ratio (Vr)/(v%) is approximately 0.66, which is a typical value in the literature. 

Next, we analyze the data with respect to scaling properties. Fig. (2^.) shows the energy spectrum with Kolmogorov's 
-5/3-law for comparison. In Fig. [3J>) the third order structure function is plotted in a compensated representation, 
i.e. (\u^\)/r is plotted against r to estimate the scaling range. The maximum lies at 10~ 2 m, which defines according 
to Kolmogorov's 4/5-law the position and the width of the scaling range. 

Because the scaling range is narrow, we use the extended self similarity (ESS) 0,0, see also appendix IVIII 

(KH cx (\u r \ 3 f (19) 



(KP> oc (K| 3 )«- (20) 

to estimate the scaling exponents of the structure functions. In Fig. [3] we show the application of ESS to the third, 
fourth and sixth order structure functions. The third order structure function is of especially interest, because it serve 
as the reference in ESS. In all three figures it can be seen that the transverse exponent is smaller than the longitudinal 
one, £* < £ l n . This result is well accepted 0, 0, 0, EH EH ■ As a consequence, one scaling group is not enough to 
characterize the turbulent flow and the statistics is more complex than previous thought. 

To make this more clear, Fig. 0] shows the exponents £} n and £* up to order 8. For the sixth order structure 
function we get ^ = 1.76 ±0.04. Fitting Kolmo gorov's log-normal model to the values of yields the intermittency 
parameter [i =0.24. Both is in good accordance with the experimentally expected values, see for example pj|. The 
transverse exponents are clearly smaller then the longitudinal exponents for n > 3, i.e. the transverse structures are 
significantly more intermittent, see Fig. This heavily discussed point will be discussed in detail at the end of this 
article. 

An other quantity which enables to quantify intermittency is the flatness F a = (a 4 ) / ( a2 )~ j where a is u r or v r , 
respectively. For a Gaussian distribution it is F a = 3 and for an intermittent distribution it is F a > 3. As shown 
in Fig. |SJ for r < L both components are intermittent. For small r the transverse increments are considerably more 
intermittent than the longitudinal one. 

In Fig. the probability density of longitudinal and transversal velocity increments are plotted. The velocity 
distribution on each scale is normalized with the respective standard deviation, i.e. u„ :— u r /a u>r and v„ := v r /a Vjr , 
to compare only the form of the curves. In Fig. \7\ the longitudinal and transversal probability densities for two 
different length scales (r = L/5 and r = A) are shown. In both figures, the intermittent character of the statistics can 
be seen. Whereas for L — L/5 the main difference is seen only for positive increments, i.e. the main difference is the 
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FIG. 2: a) Energy spectrum (dots) with a -5/3-power law and b) compensated plot of the third order structure function. The 
plateau of (|iir| 3 ) /r defines the scaling range SR. (Data for Ra = 180.) 




<v 4 > 5 P 1.05 
<u 4 > 4,=1.28 
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<l«(r)l 3 > 




FIG. 3: The ESS representation of the longitudinal and transverse structure functions, a) The transverse third order structure 
function is plotted versus the longitudinal third order structure function. The exponent is £3 = 0.87. The line has an exponent 
of 1, which is expected from Kolmogorov's 4/5-law for the third order longitudinal structure function, b) The fourth order 
structure functions and c) the sixth order structure functions. Between the longitudinal and transverse structure functions 
clear differences can be seen. 




FIG. 4: The scaling exponents estimated with ESS. Bold squares: the scaling exponents up to order 8. The dotted line 
is the forecast of K41, the straight line includes the intermittency corrections of K62 with the fitted intermittency parameter 
fi = 0.24. Open squares: transverse exponents. 
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FIG. 5: Flatness of longitudinal and transverse increments. For scales r < L the flatness lies for both components over the 
Gaussian value of 3. Furthermore the flatness is larger for the transverse component, i.e. it is more intermittent. 




FIG. 6: Probability density functions (pdf) for a) longitudinal and b) transverse increments on three different length scales 
r = L, L/5, A « A. The pdfs' width are normalized with respect to their standard deviations and are shifted along the ordinate 
for a better representation. On the largest scale L a Gaussian distribution is fitted for comparison. Towards smaller scales the 
deviations from the Gaussian form become obvious. 
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FIG. 7: A direct comparison between the longitudinal and transverse pdfs, a) for r — L/5 and b) for r = A. Longitudinal pdfs 
are represented by black squares, transverse pdfs by white squares. The width of the pdfs are normalized with respect to their 
standard deviation. For r — L/5 the deviations are mainly given by the skewness, for r — A also differences in intermittency, 
i.e. exponential tails, occurs. 
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FIG. 8: The second order transverse structure function calculated from data (squares) and calculated from the longitudinal 
structure function via the Karman equation (line). Because the differences in the scaling range can not be seen, also the relative 
error between them are plotted (dashed line). (Ra=550) 




/ [Hz] r [m] 

FIG. 9: a) Energy spectrum (dots) with a -5/3-power law and b) the third order structure function for Ra=550. The plateau 
of (|w| 3 ) jr defines the scaling range SR. 

skewness, we find for smaller scales that the distributions for the transverse increments is for positive and negative 
increments more intermittent, see Fig. 0d). 

Next we present briefly results from the high Reynolds number data set. The Karman equation is well fulfilled, i.e. 
the data are isotropic in good approximation, see Fig. [SJ Fig. [5^,) presents the energy spectrum of the longitudinal 
increments which shows a distinct scaling range with the exponent -5/3 in accordance to Kolmogorov's theory. The 
scaling range is more pronounced for this data set as it can be seen also from the third order structure function, 
compare Figs. EDd) and [2b)- 

In this section, we have presented our data with respect to two-point statistics, i.e. regarding only the statistics 
of increments for one fixed scale separately. The results are comparable to them given in the literature for moderate 
Reynolds numbers. Differences between longitudinal and transverse increments are clearly visible in pdfs and structure 
functions indicating that the transverse increments are more intermittent as the longitudinal. In the following we apply 
a new analysis to study the dependence between different length scale and the interaction between both increments 
which was not studied so far. 



B. Multi-Point Statistics: FP-Analysis 

In this section we present the analysis of the multi-point statistics separately for the longitudinal and transverse 
increments as it was described in section^] There are in essential three steps: First we show the validity of the Markov 
properties. Secondly, we calculate the Kramers-Moyal coefficients and show that the data obey a diffusion process. 
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FIG. 10: a) Contour plot of the single (straight line) and double conditioned probability distribution (dashed line) p(ui , r\ \u% , T"2) 
and p(ui, ri \u2, ri\ uz = 0,7-3) of the longitudinal increments for the length scales ri = Ar, T2 = 2Ar and rz = 3Ar with 
Ar = 68.3 mm. The distance between the contour lines is Ap = 0.05. b) Three cuts through the contour plot are presented for 
U2 = —2(7, U2 — and 112 = —2a. It can be seen that the Markov properties are fulfilled. 




FIG. 11: The same plot as in Fig. 1101 but for smaller values Ar = 2.54 mm ~ A/3, a) Contour plot for the single conditioned 
(straight line) and double conditioned probability distribution (dashed line) . Deviations between both are visible. The distance 
between the lines correspond to a factor 10. b) Three cuts through the contour plot are presented for U2 = —2a, 112 — and 
U2 — —2a. It can be seen that the Markov properties are not fulfilled. 

At last, as a verification, we integrate the resulting Fokker-Planck equation for the simple and for the conditioned 
probability distribution and show that the estimated Fokker-Planck equation describes correctly the data. 

Inserting a comment on the increment definition: different from the usual increment definition, we define increments 
for the multi-point examinations according to u r := e • [U(x + r/2) — U(x — r/2)]. In appendix IVIIII we compare 
both definitions. 



1. Markov Properties 

The foundation of the following analysis is the validity of the Markov properties. They can be tested directly 
on data by their definition J3J. Because of the finite numbers of measured data points we restrict ourself to the 
verification of p{u\, ri|«2, r^) = p(ui, 7T|it2, 7*2; U3, ra). Fig. 1101 shows both side of this equation for the longitudinal 
increments for the three length scales 7*1 = Ar, r-i = 2Ar and = 3Ar with Ar = 68.3 mmw L/2 . It can be seen 
that both distributions coincide. For a length scale below a certain threshold, the Markov properties are not fulfilled, 
see Fig. ^2 We call the associated length scale 'Markov coherence length' or short 'Markov length' l m , i.e. above l m 
the Markovian properties are fulfilled, below l m not [2^| . 

To quantify the results and to get a more objective and systematic measure for the Markov properties and the Markov 
length, we perform a Wilcoxon test, which compares two random samples with size m and n (see |54l l67j ). For the 
Wilcoxon test, one has to count the number of inversions of two samples, here for the single and double conditioned 
variable Ui| U2 and ui| U2 ,u 3 - We calculate (t) := \Q — (Q)p= p \/cr(m, n), where Q is the number of inversions calculated 
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FIG. 12: The expectation value (i) of the Wilcoxon test in dependence of the length scale differences Ar for the longitudinal 
increments with the reference scale ri = 14.64 mm; a) for M3 = and b) for 113 — a. The constant line marks the expectation 
value for fulfilled Markov properties, the dotted line is an exponential decay. 



-1 1 1 1 r- 




12 



Ar [mm] 



Ar [mm] 



FIG. 13: The same plot as in Fig. 1121 but for the transverse increments. The only differences can be found in the smaller 
Markov- length . 



from the experimental data for the variables Ui\ U2 and lfi|« 2 ,« 3 ; (Q)p= P = mn/2 and a(m, n) = ^Jmn(m + n + 1)/12 
are the number of inversions and the standard deviation, respectively, assuming that both variables have the same 
distribution. Thus, it is (t) = 1 if both samples come from the same universe, or have the same distribution. Fig. IT51 
shows this measurement in dependence of Ar. For small Ar the Markov properties are not fulfilled, whereas for larger 
distances the deviations are not significant anymore. We identify the distance Ar = l rn , where (t) drops to 1 as the 
Markov length 70] . This value can be estimated by fitting an exponential function to the values and by interpreting 
the passage through the value 1 as the Markov length, as it is presented in Fig. The Markov properties are also 
fulfilled for transverse increments, but with a smaller Markov length, see Fig. The Markov length varies within 
20% with respect to the condition ui but remains about constant with respect to r. For the longitudinal increments 
the Markov length lies in the range 7.4 mm < l m> i < 9.6 mm, for the transversal increments the Markov length lies in 
the range 5.5 mm < l m j < 6.8mm. The ratio is l m ,i/lm,t ~ 1-4 as it is known for the Taylor length [EJ. Note, that 
up to now wee see that l m ss A [54|. The analysis of the data RA=550 give in principle analogous results. 

We conclude from the results that the 'cascade' of longitudinal and transverse increments can be described by a 
Markov-process for stepsizes larger than the Markov length l m , which becomes important again for the estimation of 
the Kramers-Moyal coefficients, as will be seen below. 



2. Kramers-Moyal Coefficients 

The drift coefficients and the diffusion coefficients are calculated according to Eq. (7J directly from 

the measured data following the procedure described in [22I IH^ . The crucial point is the estimation of the limit 
limAr^o AfW, see Fig. ^1 Only the points with Ar > Z m are used to estimate the limit. For Ar < Z m the Markov- 
properties are violated, for Ar > l m the dependence of M^> with respect to Ar is linear, see Fig. ^] We thus use a 
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FIG. 14: The coefficients M' 1 ^ and in dependence of the step-size Ar for r = L/2. The two upper figures show the limit 

of the drift coefficient for the longitudinal and transverse increments, respectively. The two lower figures show the limit of 
the diffusion coefficient for the longitudinal and transverse increments. The dotted lines represent the Markov length. Circles: 
u, v = —2a, diamonts: it, V = 0, squares: u, v = 2a. 



first order polynomial to extrapolate the limit Ar — > 0. The linear dependence is the first order approximation of the 
limit, see |22j . 

The resulting drift coefficients ZK 1 ' and diffusion coefficients are shown in Fig. Elf° r the length scale r = L/2. 
The drift coefficient can be approximated by a straight line with negative slope. Small deviations from this behavior 
are visible for the transverse component. For the diffusion coefficient qualitative differences between both increments 
are visible. In contrast to the transverse coefficient, the longitudinal coefficient is not symmetric under reflection 
u — > —u. This is compatible with Kolmogorov's 4/5-law, which states that the longitudinal distribution are skewed. 
The diffusion coefficient can be approximated by a second order polynomial, so we have 

D (1) (a,r) = d?(r)a (21) 
D (2) (a,r) = d 2 {r) + d%{r)a + d% a (r)a 2 , 



where a = u, v. Due to the reflection symmetry v — » —v of the transverse increments it is d^ir) = 0. 

The drift and diffusion coefficients are the first two coefficients of the Kramers-Moyal expansion. According to the 
Pawula-theorem all higher coefficients vanish, if the fourth-order coefficients are zero and thus the expansion simplifies 
to a Fokker-Planck equation. In Fig. ^]this fourth order coefficients are plotted for r = L/2. The coefficient 
for the longitudinal coefficient vanish within the error bars. The corresponding transversal coefficient has a value 
slightly above zero. But one can estimate with the Kramers-Moyal expansion that the contribution of this coefficient 
is less than one- hundredth in comparison to the diffusion coefficient. Therefore, we want to assume in the following a 
vanishing fourth order coefficient. This assumption is additionally justified below by showing that the Fokker-Planck 
equation describes the increment statistics well. 

Next, the r-dependence in Eq. (|21|l is investigated. It can be estimated by fitting the approximation l|21|l to the 
numerical Kramers-Moyal coefficients. The result is depicted in Fig. El f° r the longitudinal (black squares) and 
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FIG. 15: The Kramers-Moyal coefficients D*- 1 ' and D^ 2 ' and their dependence of the increments u and v, respectively for 
r = L/2. a) and b) longitudinal, c) and d) transverse increments. 




u 



FIG. 16: The fourth order Kramers-Moyal coefficients D^'(u,r) and D^\v,r) for r = L/2. It can be seen that within the 
error-bars the longitudinal coefficient is zero. The transversal coefficient is slightly above zero, but its contribution to the 
Kramers-Moyal expansion is small. 



transverse increments (white squares). The form is remarkably simple, it can be approximated by 



: 


dh(r) 




: 


dUr) 




Dp) : 


d2,i(r) 


= fo,ir, dl l (r)=p% l r, dfi{r) = aft 


A (8) : 


d2,t(r) 


= for, dZ (r) = all + + ^ 



Here, we denote by X\ and X t the longitudinal and transverse quantities, respectively. 
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FIG. 17: The r-dependence of the expansion coefficients of the Kramers-Moyal coefficients according to Eq. [12 1 1 for the 
longitudinal (black squares; a = u) and transverse increments (white squares; a = v) If the notation is unique, we omit the 
index I or t. 



3. Integration of the Fokker-Planck Equation 



Next we want to demonstrate that the Fokker-Planck equation can describe correctly the statistics of the turbulent 
field. If the estimated drift and diffusion coefficients (|22[1 are used to solve the Fokker-Planck equation numerically, 
the resulting distributions can be compared with the distributions estimated directly from the data. First, the 
numerical estimation of the single distributions p(u, r) and p(v, r) in dependence of r are discussed, for which the 
distribution on the integral scale p(u,r = L) and p(v,r = L) are used as initial conditions. In Fig. 1181 v(u. r) and 
p(v,r) are shown for several length scales. The curves are in good agreement with the data. It is important to stress 
that also the intermittency effects and the skewness can be described well. A similar calculation has been done for 
the conditional distributions p(u,r\uQ,ro) and p(v,r\vQ,r ) starting at the integral scale with Dirac's delta function 
p(u,r = L\uo,ro — L) — 5(u) and p(v,r = L\vo,tq = L) = S(v) as initial condition. The solutions down to r = L/2 
are shown in Figs. 1191 and I2UI in dependence of the initial value uq and vq. Comparing the conditional probability 
distributions in Fig. ^5] with those in Fig. [5U|it is evident that the transverse statistics relaxes faster (the contour 
lines are more horizontal). 

On the basis of these results we conclude that the increment statistics can be well described by a Fokker-Planck 
equation, whose drift and diffusion coefficients are given by (|22|) . Thus, we can examine the increment statistics by 
means of the drift coefficient and diffusion coefficient D^ 2 \ 

A closer look at the d-coefficients of Eq. I|22|) shown in Fig. Elgi yes insight into the differences of the longitudinal 
and transverse increment statistic. It can be seen in Fig. 117b ) b) and d) that the absolute value of the transverse 
coefficients are larger than the longitudinal one. This means that the transverse cascade are in a way 'faster' and 
more noisy. This behavior can be state more precisely by taking into account a remarkably symmetry, namely, if the 
r-dependence of the transverse increments is rescaled by the factor 2/3, i.e. r — ► 2r/3, the dominating terms fall on 
top of the others, see Fig. [21 Thus the main statistical differences between longitudinal and transverse increments 
vanish by this rescaling. Differences only remains in the diffusion term, see Fig. [22 In section Ivl we discuss this 
result in further detail. 
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FIG. 18: The solution of the Fokker-Planck equation in comparison to the data. With the initial distribution on the scale 
r = L the Fokker-Planck equation is solved numericaly. The curves belongs to the scales r =131, 74, 41, 23, 12 mm from top to 
bottom. For a better visibility the curves are shifted by a constant factor, a) The solution for longitudinal and b) for transverse 
increments. In both pictures the intermittency is visible, for the longitudinal increments also the skewness can be seen. 




FIG. 19: The longitudinal conditional pdfs p(u2,r = L/2\u\,r = L) calculated from r — L down to r = L/2. a) The contour 
plot for the numerical solution (dashed line) and pdfs calculated from the data (straight line) . The distance between the contour 
lines is Ap = 0.05. b) Cuts through the contour plot for U2 = — 2.5<r, 0, 2.5<r. The lines are the solution of the Fokker-Planck 
equation, the symbols represent the pdfs from the data. 

C. Joint multipoint-statistics of longitudinal and transverse increments 

In the preceding section we have analyzed the statistics of longitudinal and transverse increments separately. But 
this separation is restrictive because the dynamics of both components come from one velocity field and are therefore 
connected. They can not be separated due to the nonlinear advection term in the Navier-Stokes equation. Therefore 
we extend the above analysis and examine the joint stochastic properties. We chose the both increments as one state 
with the scale parameter r, see Eq. ©. The aim is to estimate a Fokker-Planck equation © in these two variables. 
We proceed similar to the one dimension analysis. The main difference is that one need much more data and one has 
to estimate much more coefficients, because the drift-coefficient is now a vector and the diffusion coefficient a matrix. 

The verification of the Markov properties for two variables is doubtful because one has to estimate the double 
conditioned probability function for a two dimensional process, i.e. a six-dimensional function with a finite number 
of data points. One needs approximately 10 times more data points in comparison to the one dimensional case if 
the results should be similar significant. The limiting factor is the duration of the measurement and the amount 
of data. But we know separately for both components that the the Markov properties are valid. In general, if two 
variables have Markov properties, also the joint statistics have Markov properties (the reverse is not true in general 
|H7j). Therefore we assume that the combined process is Markovian as well. 

The next step is to estimate the Kramers-Moyal coefficients. First, one has to calculate the approximation of the 
drift vector m\ X '{u, v, r, Ar) and of the diffusion matrix M?' («, w, r, Ar) in dependence of Ar, see Eq. (fHJ . As in 
the one dimensional case we calculate the approximation of the drift and diffusion coefficients by fitting a linear 
polynomial to these functions in dependence of Ar above the Markov length. In doing so, we have to consider two 
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FIG. 20: The same as Fig. 1191 but for the transverse increments, a) The contour plot for the numerical solution (dashed line) 
and pdfs calculated from the data (straight line). The distance between the contour lines is Ap = 0.05. b) Cuts through the 
contour plot for U2 = 0. The line is the solution of the Fokker-Planck equation, the symbols represent the pdf from the data. 
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FIG. 21: The drift coefficient a) and the constant term of diffusion coefficient b) of the longitudinal (black squares) and 
transverse component (interrupted line) can be shifted on top of the other by stretching the r-dependence of the transverse 
coefficients by a factor 2/3 (white squares). 




FIG. 22: The diffusion coefficient of the longitudinal increments onr = L/2 (black squares) and the transverse increments on 
r — 2/3 ■ L/2 (white squares). It can be seen the minima for both are the same. 
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FIG. 23: The u and v dependence of the drift vector for the scale r = L/4. a) The drift coefficient d[ , b) the drift coefficient 
I?2 ■ Note that the vertical axis is rotated for an better comparison between a) and b). Both coefficients are pure linear and 
depend only on one variable. 




FIG. 24: The u and u-dependence of the diffusion matrix for the scale r — L/4. a) The coefficient D xl , b) the coefficient D\ 2 . 
It can be seen that the diagonal coefficient are not constant but have a parabolic form, which is more pronounced for the 
coefficient (multiplicative noise). Both coefficients are symmetric under reflection with respect to v — » — v, but not for u — > —it. 
c) The saddle- formed off-diagonal coefficient. 



different Markov length for both components, one for each component. 

The shape of the drift and diffusion coeffients are shown in Figs. [23 and [2] for the length scale L/4. The 
drift coefficients have a simple form: they depend only linearly on one variable, see Fig. 1231 That means the process 
decouples in the deterministic part. Remark if this was not the case the increments could not be Markovian seperately 
in contradiction to our findings. The diagonal coefficients of the diffusion matrix are shown in Figs. 124b ) and 124b). 
Roughly spoken they consist of a curved surface shifted upwards with a minimum at u > and v = 0. The off-diagonal 
coefficient D 12 (it is D 12 = D 21 ) has a saddle shape, see Fig. 1241). 

To summarize, we can approximate the form of the drift and diffusion coefficients by a low-order polynomial in u 
and v: 

D[ 1] (u,v,r) = <(r)u (23) 
D^(u,v,r) = <f 2 (r)v 

D^(u,v,r) = dni^+d^u + d^i^u 2 +d v 1 v 1 (r)v 2 
Dg\u, v, r) = d 22 (r) + d 2 l 2 (r)u + d^{r)u 2 + d» 2 l{r)v 2 
D${u,v,r) = d 12 (r)+d v 12 (r)v + d^(r)uv. 

The d-coefficients contain the r-dependence. Their lower index labels the associated Kramers-Moyal coefficient, the 
upper index the order of the coefficient with respect to u and v. In the Fokker-Planck equation, the coefficients 
occur symmetric with respect to reflection v — > —v. Thus a reflection v — * — v does not change anything, whereas 
this symmetry is violated for the longitudinal increment. Of course, the Kramers-Moyal coefficients can be better 
approximated using higher order in u and v, but their contributions are small and their value are not well defined 
at the edges of the available data range, because high velocities are too rare to ensure a good statistics. Note, the 
significance of higher order terms in and is important for the closure of Eqn. I|16[) but their investigation 
should not be addressed in this article. 
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FIG. 25: The coefficients d of Eq. 112MB . which reflect the dependence of the Kramers- Moyal equation on u and v in dependence 
of the linear scale variable r. In Figs, a), b) and c), the argument of the transversal coefficient is rescaled by the factor 2/3. 
The corresponding longitudinal and transversal coefficients coincides which each other apart from the intermittency term d^s;- 
e) and f): the coefficients for the off-diagonal diffusion coefficients. 



In order to describe the statistics with the Fokker-Planck equation completely, the r-dependence of the ^-coefficients 
has to be estimated, see Fig. It can be approximated by 
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Before we interpret the results, it has to be shown that the Fokker-Planck equation together with the coefficients Ij23(l 
and 124|) can reproduce the statistics of the measured data. There are in principal two ways to verify the estimated 
drift and diffusion coefficients. One way is to solve the Fokker-Planck equation, the other way is to calculate the 
structure functions, which can be done for example with the hierarchical equation for the structure functions JT^Jl. In 
both cases the results can be compared with the corresponding quantities estimated directly from the data. In Fig. 
1261 the solution of the hierarchical structure function equation Ijltjfl is shown for n = and m = 1, . . . , 6. It is in good 
agreement with the structure functions. To show that also the joint probability distributions can be reproduced, we 
solve the Fokker-Planck equation by calculating the path-integral (|10[1 numerically. We use the estimated Kramers- 
Moyal coefficients and start the simulation on the integral scale r = L with a Gaussian distribution for p(uo,r = L) 
and integrate down to r = 2A. Figs. 1771 and 1281 show the results. Because of the good accordance we assume that the 
d-coefficients can be used to characterize the statistics of longitudinal and transverse increments. 

The number of coefficients in l|24|) can be reduced by a remarkably symmetry. If one multiplies the r-dependence for 
the transverse increments with a factor 2/3, the related coefficients of longitudinal and transverse increments coincides, 
i.e. diong(^) ~ ^trans (2r/3). The only exception is the diffusion term d 22 (r). Thus the number of independent 
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FIG. 26: The even longitudinal structure functions up to order 6 calculated directly from data (symbols) compared to that of 
the Fokker-Planck equation using Eq. 1161 (solid line). 
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FIG. 27: Solution of the Fokker-Planck equation, a) Contour plot of the initial condition in logarithmic scale. The simulation 
starts at the integral length r = L with a Gaussian distribution fitted to the data, b) The contour plots in logarithm scale of 
the simulated probability distribution on the scale r = 2A. The distance between the contour lines is chosen in logarithmic scale 
and correspond to a factor 10. Dashed lines are the probability calculated directly from data, the full lines are the simulation 
ones. The simulation reproduces well the properties of the data. 




FIG. 28: Cuts through the probability distributions shown in Fig. 1271 Black circles are the pdfs of the data, the straight line 
is the simulation and the dashed line is the initial condition at r = L. a) The distribution in dependence of u for v = 0. b) 
The distribution in dependence of v for u = 0. Notice that the skewness in Fig. a) and the intermittency for both components 
are well reproduced. 
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coefficients is halved. This symmetry will be examined in more detail below. 



V. INTERPRETATION OF THE RESULTS 



From the above analysis we know that the drift and diffusion coefficients contain the information of the small- 
scale statistics of the turbulence. More specifically, they contain the joint statistics of longitudinal and transverse 
increments as well as the statistics on multi scales. Thus we can study the structure of small-scale turbulence with 
aid of the drift and diffusion coefficients. 

First we will argue that the hierarchical equation 1)16(1 is a kind of generalized Karman equation and we discuss the 
similarities and differences between these two equations. Then we consider the 2/3-rescaling symmetry in more detail, 
show that it is consistent with known results and find that the coupling of longitudinal and transversal increments 
is important for intermittency. At least we show that the scaling exponents of the transverse structure function can 
fake a too small value if the frequently used extended self similarity (ESS) is used to estimate the transverse scaling 
exponent without taking into account the rescaling. This is an important result belonging to a frequently discussed 
issue in the analysis of fully-developed turbulence: the problem of possible differences in the scaling properties of 
longitudinal and transverse velocity increments in isotropic small scale turbulence. 



A. Generalized Karman equation 

The hierarchical equation l(16|) calculated from the Fokker-Planck equation is a kind of generalized Karman equation: 
it is an equation which connects different moments. Let us compare the Karman equation 



r 



d(ul) 



dr 

with the case m = 2, n — of the hierarchy i|16|) 

d(u 2 ) 



2(u 2 )-2(v 2 ) (26) 



dr 



= (2d? + dg) (u 2 r ) + d n + dfl (v 2 r ) . (27) 



Although the Karman equation holds for Re— > oo and our result is obtained for quite moderate Reynolds number, the 
structure of both equations is remarkably similar, except for the additive term dn. In Fig. 129b ) is shown how well 
these equations reproduce — d r (u 2 )- The agreement of <|27H to the data is better than that of the Karman equation 
because the flow seems to be slightly anisotropic. Thus we can also analyze anisotropic effects with the Fokker-Planck 
equation. We want to look in more detail at the differences by comparing the terms associated to the same structure 
functions on the right hand sides of the two equations, see Fig. 129b). It can be seen that the pre- factors of the 
structure functions are different for equation l|26l) and l|27|): for example it is (2d" + d"") (w 2 ) ^ —2 (u 2 ). This 
reflects the different meaning of the two equations. Note that the anisotropy is clearer visible in Fig. 129b ) than in 
Fig. n because in 129b) the difference of two quantities of similar magnitude, right hand side of Eq. 126|) . is plotted 
and thereby the deviation to the isotropic case is more pronounced. 

Whereas the Karman equation connects only the second order structure functions, the Fokker-Planck equation 
is an equation also for higher orders. Furthermore its solution not only gives the relation between the structure 
functions but can even reproduce the structure functions itself. These informations are included in the d-cocmcients. 
An additional difference is that the Fokker-Planck equation also includes the information of the multi-scale statistics. 

We want to mention that also the second Karman equation [44| 

= <«?>- 6 (u r vl) (28) 
has a similar structure as the hierarchical equation of the Fokker-Planck equation (Eq. I|27[) for m — 3 and n = 0): 

d(u 3 ) 

-r-^r- = (3d? + 3d™) («?) + 3d^ (u 2 ) + 3d?? (u r v 2 ) . (29) 



B. Different 'cascade speed' 



In Fig. [23] we have recognized by the drift and diffusion coefficients that the dominating differences between 
longitudinal and transverse increments can be found in a 3/2-times faster 'cascade speed' of the transverse increments, 
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FIG. 29: a) Comparison of the left hand sides of the Karman equation i1'2(it with the structure function equation l|27fl derived 
from the Fokker-Planck equation. Both equations are divided by r. Open squares: — d r (m^) calculated with Eq. 1271 : crosses: 
directly from data; dotted line: calculated from the Karman equation, b) The different terms of (12711 (symbols) and I26H (lines) 
which contribute to a). The pre- factors of the structure functions of equation 1271 and 12611 are different. 



which we have interpreted as a rescaling symmetry. This has some consequences for the increment statistics. First we 
want to examine how this symmetry becomes apparent in the structure functions. Then we discuss small differences 
between longitudinal and transverse increments, which exist beyond the differences in the different cascade speeds. 
At last we show that the rescaling is compatible to known results supporting our findings. 

The rescaling-symmetry of the cascades speeds can also be observed directly at the structure functions (u™) and 
{v™). This can be studied with aid of the hierarchical equation l|16f) and the coefficients in <|23l) . We apply the 
rescaling to the transverse hierarchical equation and label all functions with a tilde, whose argument is multiplied by 
a factor 2/3. Then the equations for longitudinal and transverse structure functions read 

«"> + =*^dS «- V) (30) 

9 p< ~r^-r m(m — 1) ~ /"""ITr^ m(m — 1) ~ , ""~~^__" 2 \ 



+ ^^4>FV.). (3i) 



Due to the 2/3-rescaling symmetry the coefficients d 2 , c?22, d^u and d™ can be replaced by the corresponding coefficients 
of the longitudinal equation; the only exception is the double underlined coefficient: 

+ ^^^f) + ^=^^<^>. (32) 

Equation (|30|) and (|32|l would have the same solution without the underlined terms as can be seen by comparison. 
In other words, without these terms the longitudinal and transverse structure functions also obey the 2/3-rescaling 
symmetry, (|i> m (r)|) = (|w m (3r/2)|), were we set u r = u{r) and v r = v(r). 

To focus on the deviations of the 2/3-rescaling symmetry, we subtract l|32|l from equation l|30l) . i.e. we examine the 
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differences of (u m (r)) and (v m (2r/3)) = (v m (r)): 
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This equation is a differential equation for the differences (u™) — (u™) of the longitudinal and transverse structure 
functions. If the initial condition of this equation is zero, the solution would be zero without the underlined terms. 
Furthermore, because md\ + m(m — l)d""/2 < for all order m of interest, this equation has stable solutions and 
deviations in the initial conditions will decrease fast converging to the 2/3-rescaling apart from the discussed terms. 
The simply underlined term violates the 2/3-symmetry because of d\^(r) ^ d 22 (2r/3), the double underlined term 
violates it because of the not symmetric occurrence of the coefficient d% x in the equation. For even m the last term is 
the smallest one because of the odd moment and because of the small pre- factor, which is on the same magnitude as 
the quadratic term. The influence of the two last terms become larger with increasing order of the moments. 

Although the deviations of the 2/3-rescaling are small, they have a crucial meaning for the small scale turbulence: 
without the terms c?"" and d 22 there were no intermittency, without the term d^ there was no skewness of the 
longitudinal increments. To study how the form of the Kramers-Moyal coefficients leads to increasing intermittency 
with decreasing scales, we use the flatness F u = (uf) / (it 2 .) . For a Gaussian distributed function, the flatness has 
the value 3 and deviations from this value can be interpreted as intermittency. If we differentiate the flatness with 
respect to r, we can relate the Fokker-Planck equation and the flatness (we shorten S mn := (u™v n )): 

_™11 [d c20c02 qc22] o 1 I"oc30,jtt , qg3Qju 1 
(£20)2 L^"* 5 ^ 6A J (520)2 l ZA «l+' 3 < : ' "llj 

r^- = -4F v d™ + ^ [F v - 3] + (35) 
F V S 02 S 20 - 3S 22 ] - 2—^ [2S 03 d v 2 + 3S 21 d% 2 ] . 



(^02)2 L v (502)2 

The similar structure of Eq. i|34[l and (|36|l allows to discuss both equations together. The direction of the cascade is 
from large to smaller scales r, so that d r F u < means increasing intermittency towards smaller scales. The first term 
is the dominating negative term. The second term vanishes for large scales, because of the approximately Gaussian 
form of the distribution and is positive if intermittency increase. The third term simplifies to j^xrjz0~ 2 2 because of the 

approximate Gaussian character for large scales with the covariance a\ 2 . This term describes the coupling between 
both components and is larger than zero, i.e. it acts against intermittency. The fourth term belongs to higher order 
corrections and describes the influence of the skewness onto the intermittency. To summarize, only the two quadratic 
terms dff and d 22 significantly produce intermittency of the longitudinal and the transverse component, respectively. 

Similar consideration can be done for the skewness S u := (uf) /((it 2 .)) 3 / 2 , which can also be expressed by the 
Kramers-Moyal coefficients (it is S v = because of reflection symmetry) : 

dS u 3 S 3 3 f 1 S 02 S 12 \ 

r ~d7 = 2 5^ dn ~ (520)1/2 d "i ~ 2 Sud i" + 3Su [2 S 2 ^ ~ J ^ 

A positive derivative d r S u > means an increasing skewness for decreasing scales. The first term on the right hand 
side of Eq. I|3()l) is always negative. The second term is positive and is therefore necessary for a skewed distribution. 
The third term amplifies an existing skewness. Therefore the coefficient d^ is essential for the skewness, but also the 
'intermittency' term d%± leads to an increasing skewness. 

We can summarize the results in a short form: 
If the r -dependence of the transverse component is reseated by a factor of 2/3, then the only differences can be found 
in intermittency contributions and the skewness. 
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FIG. 30: a) The transverse second structure function calculated on the one hand by the Karman equation (straight line) 
and on the other hand by using the 3/2-rescaling (broken line). Both curves deviate only within 4% from another up to 
the integral scale. For comparison the longitudinal structure function is shown (dashed line), b) The local exponent for the 
transverse structure function calculated from the longitudinal structure function using the Karman equation (straight line) and 
the 3/2-rescaling (broken line). The exponents are almost the same. 



C. Compatibility to known results 

To support our findings of the rescaling symmetry, we show that it is compatible to the Karman equation and to 
the ratio of Kolmogorov constants. These two aspects were already presented in |60| . but here we go in more detail 
and add some new aspects. 

The Karman equation 1)18(1 is consistent with the 2/3-rescaling symmetry, as one can see by interpreting the Karman 
equation as a first-order Taylor expansion with the 'small' quantity r/2: 

{v 2 (r)) = ( u 2 {r)) + l -r^-{u 2 {r)) (37) 

u\r+\r)^ - R= (u\lr)^ - R. 

This equation correspond to the 2 /3-rescaling symmetry of the second order structure functions except for the La- 
grange remainder [ljj R — r 2 d pp (u 2 (p)) /8 with r < p < 3r/2. If we assume Kolmogorov-scaling |37| (u 2 ) oc r 2 / 3 
(intermittency effects for the second order structure function are negligible), the remainder R — r 2 ^(C2 — l)p^ 2 ~ 2 /8 
is below 2.1% of (u 2 (3r/2)) for all r. An exact calculation yields that the relative error R/ (u 2 (3r/2)) is independent 
of r and R is w 1.7% of <'u 2 (3r/2)), a value below the typical errors for (u 2 ) of data from hot-wire anemometry. For 
the used data the approximation of (u 2 {r)) + rd r (u 2 (r)) jr by (u 2 (3r/2)) is valid within 4% for all scales, see Fig. 
130b ). Thus the 2/3-rescaling is a very good approximation of the Karman equation. Because the only assumptions for 
these equation is a solenoidal and isotropic field, the 2/3-rescaling is a property of isotropic turbulence, which should 
hold also for the limit Re^ oo. 

It is important to discuss the quality of the approximation (|37|l with respect to scaling exponents. The remainder 
R = r 2 /8^2(?2 — l)/^ 2-2 vanishes identically only for structure functions linear in r (£2 = !)■ But as one can see from 
(u 2 (r)) +r/2 9 r (« 2 (r)) oc r^ 2 and (u 2 (3r/2)) cx r^ 2 the exponent is neither changed by the above approximation 
for pure scaling nor is the local exponent changed significantly for real data as one can see from figure ISUb). The 
reason is that the validity of the approximation depend only weak on the exponent. In reverse this does mean that 
the 2/3-rescaling does not distort the exponent strong. 

The 2/3-rescaling can be related to the ratio of the Kolmogorov constants. Let us suppose that the structure 
functions scale with a power law, (v n (r)) = c"r^" and (u n (r)) = c™r^™, even though our measured structure functions 
are still far away from showing an ideal scaling behavior |5q. Note that the c" constants are related to the Kol- 
mogorov constants but includes the energy dissipation ^e"^ 3 ^> and the enstrophy (£!™^ 3 )- If we neglect the differences 
in intermittency and skewness, we can relate the structure functions according to the above mentioned rescaling: 

(v n (r)) = (u"(3r/2)) = c?r& = cf (3r/2)*». We end up with the relation £} n = & and c"/cj l = (3/2) C ". For n = 2 
and n = 4 we obtain cf/cf w 1.33 and cf/cf ~ 1.72, which deviates less than 3% from the value of c 2 /cf = 4/3 



24 



A 1.4- 



■^1 1 1 I I I I I | 1 1 1 I I I I I | 



v 1.0- 



0.6 1 1 — i — 

10' 3 10" 2 



_l I 



10" 



2.0 
1.5 

s 
V 

^ 1.0 

a" 

v. 0.5|- 
0.0 



_ b) 

■^tttm ™"""nnnnnnri_ E 


1 1 

□ _ 
□ - 

□ 

□ 

□ 

□ 





integral 

scale 
i 



r [m] 



0.0 0.1 0.2 0.3 0.4 0.5 

<m 2 >, <v 2 > [m 2 /s 2 ] 



FIG. 31: a) The ratio (v 2 ) / (v 2 ) (squares) calculated from data and the expected ratio if the structure functions would show 
Kolmogorov-scaling (straight line). The dashed line shows the isotropic case (^ 2 ) iso / ( - " 2 )i where (^ 2 ) iso is calculated with the 
Karman equation. The just small differences between the two curves indicate that the deviations from the straight line can not 
be explained with anisotropy. b) The almost constant ratio r((u 2 ) = s)/r'((u 2 ) = x) of the 'inverse' of the structure functions 
close to the value 2/3 indicates that the differences between both structure functions can be found in their scale dependence 
(see text). 



and cf/cf = 16/9 given in pj. Thus one needs only one constant 2/3 to explain the two ratios cf/cf and cf/cf. We 
conclude that the 2/3-rescaling is the underlying relation between longitudinal and transverse Kolmogorov constants 
and therefore gives a forecast also for Kolmogorov-constants of higher orders. 

For finite Reynolds number the differences between the ratio of the Kolmogorov constants and the 2/3-rescaling 
become more obvious. The ratio (w 2 ) / (u 2 ) — (c"r^») / '(cfr^) = Ct/c\ = 4/3 is not fulfilled anymore, see Fig. 
131b ). Whereas for high Reynolds numbers this quantity can be interpreted as a ratio of amplitudes of the structure 
functions, this meaning gets lost for finite Reynold numbers. But if we treat the structure function as the independent 
variable and r as the dependent variable, i.e. r is a function of the structure functions, and calculate the ratio 
r((w 2 ) = x)/r'{(u 2 ) — x), we get an almost constant value close to 2/3, see Fig. 15Tb). This is a remarkably property: 
The constant ratio of the Kolmogorov constants are just a special case of the 2/3-rescaling. Whereas the former 
is valid only for high Reynolds number, the 2/3-rescaling is also a property of moderate Reynolds- number flows. A 
further consequence belongs to the discussion whether the Kolmogorov constants are universal or not |6J, |52], |68j . Even 
if the Kolmogorov constants are not universal, the ratio r((|w m |) = x)/r' ((\u m \) — x) seems to be universal if the 
turbulence is isotropic. 



D. Intermittency and the coupling of longitudinal and transverse increments 

From the form of the Fokker- Planck equation it can be seen that intermittency is closely related to the interaction 
of longitudinal and transverse increments. More precisely, a necessary condition of intermittency is the mutual 
dependence of longitudinal and transversal increments. This can be deduced as follows. 

First, we neglect the small differences in the intermittency and the skewness. If the 2/3-rescaling is applied, the 
form of the Fokker-Planck equation becomes symmetric with respect to exchange of u r and v r . Moreover, because 
of the linear drift term and the rotational symmetric diffusion term, the Fokker-Planck equation becomes rotational 
symmetric and one can use circular coordinates with radius p 2 = u 2 + v 2 . Thus rotational symmetric functions 
F(p 2 = u 2 + v 2 ) solve the Fokker-Planck equation. 

To consider the question wether u r and v r are coupled, we construct a contradiction by assuming that u r and 
v r are independent, in the sense that F(u 2 ,v 2 ) — F(u 2 )F(v 2 ). This means for a solution F the relation F(p 2 ) = 
F(u 2 + v 2 ) = F(u 2 )F(v 2 ) holds. But this equation is only fulfilled for Gaussian distributed functions F what is in 
contradiction to the fact that the turbulent flow is intermittent. Thus we conclude the two increments u r and v r 
depend of each other. Intermittency models has to take into account this coupling between the two increments and 
the one dimensional intermittency models has to be extended. 
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FIG. 32: The ESS representation of the longitudinal structure functions in comparison to the ESST representation of the 
transverse structure functions, a) The ESS representation of the transverse third order structure function. The exponent is 
£3 = 1.02 and with that close to Kolmogorov's 4/5-law, see the line with the exponent 1. b) The fourth order structure functions 
and c) the sixth order structure functions. The transverse structure functions come significant closer to the longitudinal ones 
in comparison to Fig. OH 



E. Transverse scaling exponents 

In this section we focus on the wide debate about possible differences between the scaling exponents of longitudinal 
and transverse structure functions and consider the implications of the rescaling symmetry on it. In Fig. OH we have 
shown the recognized result that using ESS the transverse exponents are smaller than the longitudinal. Nevertheless, 
we use a new ansatz to reconsider this result and combine the differences of the longitudinal and transverse cascade 
speeds with ESS. 

ci ..... 

Because the assumption for the longitudinal structure functions in ESS is (\u n (r)\) oc (rfi(r)) n , see (|4Up. the 

corresponding relation for the transverse increments has to be (|i> n (?")|) oc {rft{r)) with a function f t different from 
Notice that the implicit assumption in ESS was ft = fi, see (|41l) and l|42|l. Using for the second order structure function 
the 3/2-rescaling found from the Fokker- Planck equation and assuming that the differences between skewness and 
intermittency are small for this order, we find a relation for the two functions: Starting with (|w 2 (?')|) = (|u 2 (3r/2)|) 

we get rft(r) = (3r/;(3r/2)/2) 2 2 (we have chosen the proportionality in ESS without loss of generality in such a 
way, that the equals sign holds). Because the intermittency corrections is small for the second order, it is £ 2 ~ £2 anc ^ 
we get 

ft(r) = 3/,(3r/2)/2. (38) 

For arbitrary structure functions it then holds 

(\v(r)\ n ) « (|r/x(|r)J " <x (\ u ^ r )\ m f^ = (\u(^r)\ 3 )^ (ESST), (39) 

which we call in the following ESST, extended self similarity for the transverse structure function. 

Fig. 13 21 shows the application of ESST to the transverse structure functions. As an remarkably result, the differences 
between both exponents vanish. Thus one scaling group is sufficient to characterize both increments. 

Notice that the differences between the exponents found with ESS are due to a none existing scaling behavior of 
the structure functions with r. It is evident that our rescaling does not change the exponents in case of pure scaling 
behavior (|u(r)| n ) oc (3r/2)* n oc r* n , which is expected to be valid if the Reynolds number goes to infinity (see also 

m- 

Next we perform the analog analysis for our second data set with the higher Reynolds number R\ = 550. In Fig. 
133b ) the sixth order structure function is plotted using ESS and in Fig. 133b ) using ESST. Again we find that the 
differences in the scaling exponents vanish for ESST (ESS: $ = 1.74 ±0.03, £| = 1.60 ±0.03; ESST: $ = 1.74 ±0.03, 
£g = 1.75 ± 0.03). In Fig. 03 the exponent ^ as well as is plotted up to order 8. The differences between the 
exponents ESS vanish within the uncertainties applying ESST instead of ESS. 

To see in more detail the influence of ESST to the exponents, in Fig. OS] the local exponents £"(r) = 
<91og(a")/<91og(|u| 3 ) (a = u, v n = 4,6) are shown. The use of ESS for the longitudinal exponents results in 
an almost constant local exponent. The transverse exponents have oscillations, which were explained by log-periodic 
oscillations. If ESST is applied, the differences between longitudinal and transverse exponents vanish within these 
oscillations. 
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FIG. 33: The sixth order longitudinal (black squares) and transversal (withe squares) structure function in a) ESS representation 
and b) ESST representation. 
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FIG. 34: Scaling exponents up to order 8 for longitudinal and transverse structure functions. Dashed line: longitudinal 
exponents. Open squares: transverse exponents estimated with ESS. Filled squares: transverse exponents estimated with 
ESST. 
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VI. CONCLUSION 



We have analyzed the differences and similarities between longitudinal and transverse increments. Differences in 
their structure functions and probability distributions were found, which were observed already from other groups. To 
look in more detail at these differences we have used a method to estimate by pure data analysis a phenomenological 
Fokker-Planck equation and have extended this analysis to both, longitudinal and transverse increments. We have 
checked carefully the mathematical prerequisites of this method and found that they are well fulfilled for our data. 
Then we have estimated the drift and diffusion coefficient of the Fokker-Planck equation directly from the data. The 
solution of the Fokker-Planck equation with these coefficients reproduces well the probability distributions and the 
structure functions. Thus, the statistics of the joint probability of longitudinal and transverse increments is encoded 
in the drift and diffusion coefficients. These new quantities for analyzing the turbulent flow determine a stochastic 
process in r, which can be interpreted as a 'cascading process'. The statistics of both increments are decoupled in the 
deterministic drift vector but are coupled via the nonlinear and no diagonal diffusion matrix. 

A remarkably result is that in a first approximation the longitudinal and transverse drift and diffusion coefficients 
can be transformed to each other by a simple rescaling, namely, by multiplying the scale of the longitudinal increments 
with the factor 3/2. With this rescaling the frequently discussed issue of the differences between longitudinal and 
transverse structure functions can be explained. In this context the extended self similarity method (ESS) has 
to be changed for the transverse structure functions according to the rescaling symmetry, leading to an extended 
self-similarity for the transverse increments (ESST). From it, we get the remarkably result that the previous found 
differences in the scaling exponents vanish. Accordingly, the seemingly differences between the exponents were a 
result of the non-scaling behavior. Thus it is evident that with increasing Reynolds number the differences in the 
scaling exponents diminish, as more and more ideal scaling behavior occurs. It is interesting to note that the 3/2 
rescaling, which we interpret as different speed of the longitudinal and transverse cascade, is compatible with the 
Karman equation. This leads to the proposal that the 3/2 rescaling is valid for large Reynolds number Re-^ oo. 

Besides the features of structure functions our method of analyzing the statistics by means of a Fokker-Planck 
equation, provides further insights in the complexity of the cascading process. Access to stochastic properties on 
multi-scales and the interaction between both increments are given. In this closer look principal differences between 
longtidudinal and transverse increments can be identified, which are not grasped by the investigation with structure 
functions. The physical meaning of these differences has to be clarified for an extended understanding of the small 
scale turbulence. 

We acknowledge fruitful discussions with R. Friedrich, A. Naert and teamwork with S. Luck and F. Durst. This 
work was supported by the DFG-grant Pe 478/9. 



VII. APPENDIX A: EXTENDED SELF SIMILARITY 

Benzi et. al. 0,0 have noticed that structure functions show extended self similarity (ESS) and that this property 
can be used to estimate scaling exponents for moderate Reynolds number for which a scaling regime is not well- 
developed. The basic assumption of ESS is that the structure functions of different order have a similar shape, which 
can be written as 

(KD ex (r/(r))«» (40) 

with the common function f(r). Additionaly it was shown that (|u 3 |) and (it 3 ) have the same scaling properties and 
if one choses according to Kolmogorov 's 4/5-law C| = 1 

one can write ESS in the form 

(Kr>cx(H 3 > ct . (4i) 

Using this expression, the scaling range extends down to 3r]-5r] and the exponent is universal, i.e. does not depend 
on the Reynolds number. 

For the transverse structure functions ESS was applied in two different ways. Some groups have plotted (v™) oc 
(|v r | 3 )^» [13, H^ . But in recent time it has been argued that 

(M") oc (K| 3 )€™ (42) 

is theoretical more justified because Kolmogorov's 4/5-law gives an exact pre diction for the longitudinal third order 
structure function and therefore (|w r | 3 ) is a good point of reference [3. fl5L fill l49l l69j . 
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FIG. 36: Properties of the left-bounded increments, a) The same plot as in Fig. 1121 but for left-bounded increments. For the 
case m = a the Markov properties are violated, b) The solution of the Fokker-Planck equation (dotted line) in comparison 
to the conditional probability distribution (full line) estimated from data. The conditioned pdfs can not be described by this 
Fokker-Planck equation. 



VIII. APPENDIX B: LEFT BOUNDED INCREMENTS 



In the case of multi-point statistics with n increments on different scales ri, . . . , r n , one has to define the relative 
position of increments. A general definition is u r (a) = e • [U(x + or) — U(x — (1 — a)r)] with a S [0, 1]. We call 
the case a = 1 left-bounded increments and the case a = 1/2 centered increments. The left-bounded increments are 
the usual way to define the relative position of several increments, see for example |37) . Although, for this increment 
definition, the Markov properties are violated as one can see in Fig. 136b) . For = the Markov properties are 
fulfilled. For 113 = a we get Markov properties around the expected Markov-length, but for larger scales the Markov 
properties are violated. Nevertheless, the violation of the Markov properties does not affect the Kramers-Moyal 
coefficients very much; even the deviations in does not effect the total Fokker-Planck equation very much, because 
this term is small in comparison to the others, see Fig. 1371 The violated Markov properties become more obviously, 
if the conditioned moments estimated from the Fokker-Planck equation are compared with those from the data, see 
Fig. 136b). It can be seen that the left-bounded increments does not describe the conditional moments. Therefore we 
use the centered increments for our analysis. For a discussion on the background of random numbers and random 
walks see Ref. 66 
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the deviations are significant, see c) and d)). But because these two coefficients have only a small contribution to the drift and 
diffusion coefficients, the effective influence on the Fokker-Planck equation is nevertheless weak. 
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